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ABSTRACT 

We present time series measurements of chromospheric activity for more than 
2600 main sequence and subgiant stars on the California Planet Search (CPS) 
program with spectral types ranging from about F5V to M4V for main sequence 
stars and from GOIV to about K5IV for subgiants. The large data set of more than 
44,000 spectra allows us to identify an empirical baseline floor for chromospheric 
activity as a function of color and height above the main sequence. We define AS 
as an excess in emission in the Ca II H&K lines above the baseline activity floor 
and define radial velocity jitter as a function of AS and B — V for main sequence 
and subgiant stars. Although the jitter for any individual star can always exceed 
the baseline level, we find that K dwarfs have the lowest level of jitter. The lack 
of correlation between observed jitter and chromospheric activity in K dwarfs 
suggests that the observed jitter is dominated by instrumental or analysis errors 
and not astrophysical noise sources. Thus, given the long-term precision for 
the CPS program, radial velocities are not correlated with astrophysical noise 
for chromospherically quiet K dwarf stars, making these stars particularly well- 
suited for the highest precision Doppler surveys. Chromospherically quiet F and 
G dwarfs and subgiants exhibit higher baseline levels of astrophysical jitter than 
K dwarfs. Despite the fact that the rms in Doppler velocities is correlated with 
the mean chromospheric activity, it is rare to see one-to-one correlations between 
the individual time series activity and Doppler measurements, diminishing the 
prospects for correcting activity- induced velocity variations. 
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Introduction 



During times of high solar activity, the flux in an extended network of calcium lines 
brightens significantly in narrow band filters centered on the near-UV Ca II lines. This 
observation led O. C. Wilson to anticipate that measurements of emission in the cores of 
Ca II lines would provide a good index of stellar chromospheric activity and h e began a sur vey 
at Mt. Wilson Observatory to search for analogs to the solar sunspot cycle Jwilsonlll968l ) in 
bright nearby stars. Wilson developed a protocol for measuring Shk values, defined as the 
ratio of flux in the Ca II line cores to flux in nearby continuum regio ns. These Shk values 
are now well established indices for estimating chromospheric activity. iDuncan et al.l (119911 ) 
published seasonal Shk for 1296 stars o bserved as par t of th e Mt. Wilson Observatory 
(MWO) project from 1966 to 1983 and iBaliunas et al.l (119951 ) monitored a subset of 111 
MWO stars to search for Maunder Minimum cycles. Other la rge scale surveys h ave also 
been carried out to monitor variability in chromospheric activity. iHenry et al . (1996) derived 



Shk y alues for more than 800 nearby stars from the southern hemisphere; IStrassmeier et al. 



(200 0J) mo n itored Shk as a tracer for rotation period to find candidates for Doppler imaging; 
Hall et al.l (120071 ) have measured H and K fluxes since 1994 to search for stars with sunlike 
spot cycles. These surveys calibrate coefficients for the continuum and H & K line cores so 
that the Shk values are on a scale that matches the original MWO Shk project. 

To account for different continuum flux levels near the Ca II lines for stars of different 
spectral type, the Shk values are often parameterized as log-R HK , a logarithmi c fraction of 



the flux in the H and K line cores to photometric contributions from the star (INoyes et al. 
1984 ; iMiddelkoopI Il982[ ) . Because chromospheric activity declines with stellar age, cluster 



stars have been used to calib rate the log.R HK yalues to the rotation periods a nd ages for 
stars with 0.4 < B-V < 1.0 JNoves et al.lll984t iMamaiek & Hillenbrand! l2008h . 



Stars in the field rotate more slowly and have less chromospheric activity than their 



young cluster counterparts 
Simon. Herbig fc Boesgaard 



1985 



Wilson! Il963l; iKraftl Il967l; ISkumanich 
~ Simon fc Drakel Il989h . 



1972 



Soderblom 1983 



Gray k Nagarl (119851 ) find that 



most subgiants redward of G0IV also rotate slowly and have low chromospheric activity. 
Simon &: Drakel (119891 ) measured UV emission as an activity indicator for a set of stars with 
masses between 1.2 and 1.6 M & . They find that UV activity levels are essentially constant 
as lower mass stars with M < 1.25M evolve from the main sequence across the subgiant 
branch. In contrast, stars with masses in the range of 1.25 - 1.5 M were found to exhibit 
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moderate levels of activity (relative to the lower mass subgiants) on the subgiant branch 
with an abrupt drop in UV emission (activity) near spectra l type GOIV that was e mpirically 
accompanied by a decrease in stellar rotational velocities. ISimon &: Drake! ( 119891 ) interpret 
this activity-rotation boundary as a physical transition from acoustic heating in early F- 
type main sequence stars to magnetic dynamo-driven activity followe d by strong rota tional 
braking and a decrement in chromospheric activity. As a case in point. lGillilandl (119861 ) notes 
unusual high vsini and moderate chromospheric activity in the F8IV star, HD 125840 and 
interprets this as a stage in subgiant evolution that might develop with the onset of intense 
dynamo activity as the convective zone deepens. 

Chromospheric activity is correlated with spots and flows in stellar atmospheres. High 
chromospheric activity is empirically associated with both rando m and quasi-periodic veloci- 



2002; 


Queloz et al. 


2001; 


Bonfils et al. 


2007) 



stars that may have photospheric features that produce radial velocity "jitter." Chromo- 
spheric activity is monitored in the stellar spectra used to deriye pre c ision radial velocities 



for stars on many Doppler planet searches (ISaar. Butler. Marcv 



Paulson et all 12002 : IWright et all 12004 ; IWright 



2005; 



1998 



Johnson et al. 



Saar fc Fischerlll99&! 



2007m . IWrightl (120051 ) 



measured Suk f° r t ne California Planet Search (CPS) observations obtained before 2004 
August with the original HIRES CCD at Keck Observatory. 

Here we provide updated activity measurements with a calibration of Shk to the MWO 
scale for CPS spectra of main sequence and subgiant stars obtained after the upgrade of 
the HIRES CCD in 2004 August and activity measurements for stars on the Lick planet 
search project. We define an excess in Suk above a baseline floor (defined empirically as 
the tenth percentile level in chromospheric activity) as AS for stars for the color range 
0.4 < B — V < 1.6 and examine correlations between AS and radial velocity jitter. 



2. Keck Observations and Analysis 



The CCD detector on Keck HIRES (IVogt et all [l99J) was replaced in August 2004 
with three new detectors that span a broader wavelength range and have higher quantum 
efficiency (particularly at blue wavelengths) and improved charge transfer efficiency. The 
HIRES B5 decker is generally used for CPS program observations and provides a spectral 
resolution of about R = 52, 000. The spectral format of HIRES is configured by adjusting 
the echelle and cross- dispersor angles so that (for CPS program observations) the iodine 
spectrum falls on the middle of the three chips. At the beginning of every CPS observing 
run, the focus algorithm is used to recenter a set of Thorium and Argon emission lines so 
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that the lines always fall on the same CCD pixels. This setup procedure ensures that stellar 
absorption lines fall on nearly the same part of the detector (modulo barycentric and stellar 
velocity shifts). This provides more consistency in the line spread function, relevant for high 
precision Doppler analysis. 

The blue CCD detector spans 3650 - 4790 A and includes the Ca II H & K lines 
used to assess chromospheric activity. Program observations are made through the iodine 
absorption cell, however, the molecular iodine lines disappear blueward of 5000 A and do 
not contaminate the Ca II H&K lines. This makes it possible to serendipitously measure 
stellar activity with every Doppler velocity measurement. 



2.1. Keck Shr Calibration 



Following a method described in iDuncan et al.l (119911 ). we measure flux in the core of 



the Ca II H & K lines relative to continuum bands. Figures 1, 2 and 3 show the echelle 
orders used to measure S'hk from a spectrum of the bright G8V star r Ceti. The core flux 
in the Ca II lines is weighted using a triangle function with a FWHM of 1.09 A centered 
at 3968.47 and 3933.66 A for the H and K lines, respectively. The triangular weighting 
function gives the greatest weight to the core of the Ca II lines, which are formed higher in 
the chromosphere, and diminishing weight away from the line core. Flux is also measured 
in two continuum windows. The continuum window on the redward side of the H and K 
lines is labeled as the R section and the continuum window blueward of the lines is defined 
as the V section. The R and V continuum sections monitor the overall flux of the star 
and help account for variable observing conditions such as changing air mass. As described 
below, coefficients for flux measurements in the line core and continuum windows are free 
parameters that scale the Keck data to match Shk values from the Mt. Wilson H&K 
project for stars in common to both projects. 

Because the wings of the Ca II lines stretch across a significant fraction of the echelle 
orders, it is difficult to define the continuum near these lines. As a substitute for continuum 
normalization, we construct a template for each star using five high signal-to-noise spectra 
for each star. Each of the five stellar spectra are shifted to rest frame wavelengths by cross- 
correlation with the National Solar Observatory (NSO) solar atlas. A median filter is applied 
pixel by pixel in the stack of five observations to smooth over cosmic rays producing a high 
quality template. The time series spectra are then compared to the template, enabling 
precise measurement of variation in the Ca II line cores. 



For each of the time series spectra (program spectra, obtained for Doppler measure- 
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merits), the echelle orders with the Ca II lines are cross-correlated and shifted to the tem- 
plate wavelengths. The program spectra are then scaled in intensity to approxiately match 
the continuum intensity of the template spectra near the H & K lines. The program spectra 
generally have different continuum slopes than the template because the barycentric veloc- 
ity will differentially shift each of the time series spectral lines by several pixels. The blaze 
function of the echelle produces a maximum in flux near the center of the echelle orders that 
declines toward the red and blue edges of each order. As a result, the intensity is differen- 
tially affected by the steep blaze function across the order. To align the continuum, each 
wavelength-shifted program spectrum is divided by the template spectrum and the residual 
is fit with a second degree polynomial. When aligning the continuum, the cores of the Ca II 
lines are masked out so that any real variation in the line core flux is not affected. The 
program spectrum is then divided by the polynomial so that it is finally well-aligned with 
the template both in wavelength and continuum intensity. 

After this careful alignment with the template spectrum, small relative differences in 
the Ca II line cores can be measured to derive time series S'hk values. Because these mea- 
surements of chromospheric activity are simultaneous with the Doppler measurements, they 
provide a valuable diagnostic for interpreting prospective dynamical velocities in our Doppler 
planet survey. 

The calibration of coefficients for the continuum and core flux values was established 
using 151 main sequence stars on th e CPS Keck survey that also have MWO S'hk values 



published in the iDuncan et al.l (119911 ) catalog. The calibrating sample was restricted to stars 
with B — V < 1.5 that sit within one V magnitude of the main sequence. We initially fit 
four scaling coefficients (Eqn 1) using a Levenberg-Marquardt (LM) algorithm: 

{CiH + C 2 K) 

Srk ~ (C 3 R + C 4 V) (1) 

This produced covarient coefficients that were extremely sensitive to the initial guesses, 
so we reduced th e numbe r of fr ee co efficients by fixing the r atio of flux meaurements in the 



H and K bands. IWilsonl (119681 ) and IShine fc Linskyl ( 11974] ) both note that the flux in the 



H and K lines should be equal in intensity, and that deviations from this are instrumental. 
Examining our spectra we find an average ratio of the H line flux to the K line flux of 1.45. 
In other words, multiplying the K line flux by 1.45 boosts the K line intensity so that it is 
comparable to the H line. Similarly, we found that flux in the redder R continuum window 
was dramatically higher than flux in the blueward V continuum windows. We adopted a 
fixed coefficient of 25 to adjust flux in the V continuum window. Fixing the coefficients for 
K and V breaks the covariance between the remaining coefficients and simplifies the form of 
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our calibration to Eqn (2). The two coefficients in Eqn (2) were determined with Levenberg 
Marquardt fitting: 



S: 



_ (H+1.45K) 
SuK ~ Cl (R + 25V) + ° 2 

(H + 1A5K) 
32.510^ — z^ttt 1 + 0.021 



HK 



(R + 25V) 



(2) 
(3) 



In Figure 4 we plot Keck Shk va lues derived with Eqn (3) against standard Mt. Wilson 
Shk values from iDuncan et al.l (Il99l[ ). We find an rms scatter of 11% to this fit that is likely 
the result of both measurement uncertainty (from MWO as well as from the Keck data) and 
intrinsic chromospheric variability. Chromospheric activity of stars can vary on time scales 
of days to years. An RMS scatter of 11% means that the Shk values published here are 
within 11% of the long term average Shk of the MWO calibration stars. 

To empirically assess our measurement errors, we examined the RMS scatter for the 
chromospherically quiet, old G8V star, r Ceti. This star has a large number of observations 



(1991 


) and 


Wright et al. 


(2004) 



in Figure 5. The full width half maximum (FWHM) of this distribution is about 0.002, our 
empirical assessment of measurement error. Thus, the measured Shk at Keck for five years 
of data for r Ceti is 0.167 ± 0.002. 

Program spectra that have a signal-to-noise < 5 in the continuum near the Ca II lines 
exhibit dramatically increased scatter in the measured Shk values. We adopt a minimum 
signal-to-noise of 5 as a threshold for Shk measurements. Fewer than 2% of the stellar 
spectra analyzed here had signal-to-noise below this rejection threshold. 



3. Lick Observations and Analysis 

The Lick Observatory Planet Search began in 1987 with the Hamilton Spectrometer 
(Vogt 1987) and both the Shane 3- meter telescope and the 0.6-m Coude Auxiliary Telescope 
(CAT). The CCD detector at Lick was changed in 2002 to a Lawrence Berkeley Laboratory 
high resistivity CCD and the spectral format was extended at that time to include the Ca II 
lines for simultaneous monitoring of chromospheric activity. The current spectral format 
ranges from 3800 to 9000 A with a spectral resolution of R = 55, 000 at 6000 A . The 
quantum efficiency (QE) of the detector is about 80% for wavelengths between 5000 - 6000 
A where the iodine absorption lines are analyzed to measure Doppler shifts. The QE drops 
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to less than 50% near the blue Ca II lines. The spectra from the smaller aperture telescopes 
at Lick typically have lower signal-to-noise than Keck, particularly near the Ca II lines. 



3.1. Lick Sh Calibration 



The chromospheric activity measurements at Lick Observatory only make use of the 
Ca II H line (Figure 6) because of inadequate signal-to-noise near the echelle order containing 
the Ca II K line. Therefore, only an Sh value is cal c ulated for spectra obtained at Lick 
Observatory. We adopt the approach of IWright et al.l (120041 ) and measure the ratio of the 
line core flux in the Ca II H line relative to a single continuum window, C: L = H/C. 

Similar to our analysis of Keck stars, a coadded, median-filtered, template is created for 
each of the stars on the Lick program. One additional step is made to eliminate cosmic rays 
or pixels affected by Compton scattering of electrons in the CCD for program observations. 
This is more of a problem at the Lick Observatory because of the longer exposure times and 
a thicker substrate in the high resistivity CCD, which is more subject to Compton scatter- 
ing. To identify affected pixels, the wavelength-shifted and flux-scaled program spectrum 
is divided by the template observation. Pixels with values that are ten sigma away from 
the median value are replaced and the program spectrum is then iteratively realigned with 
the template. Spectra with replaced pixels in the continuum window or in the line core 
are flagged and excluded if they are outliers so that cosmic ray cleaning does not affect the 
Sh measurement; it simply improves the alignment of the program observations with the 
template. 

Coefficients for the L ick Sh values (Eqn 4) were calibrated using 83 stars in common 
with iDuncan et al.l (Il99ll ). Only stars with a color B — V < 1.6 that were within one V 
magnitude of the main sequence were used in the calibration. A Levenberg Marquardt fitting 
analysis was used t o find the coef fi cients that minimized the RMS scatter between the Lick 
Sh values and the IDuncan et al.l (Il99ll ) Shk measurements. This calibration has an rms 
scatter of 11%, shown in Figure 7. We tested inclusion of higher order terms but did not 
find improvement in the rms fit. 



S H = 2.206L 2 + 6.907L (4) 

To estimate our Sh measurement uncertainty at Lick Observatory, we again considered 
time series Sh measurements for the chromospherically inactive star, r Ceti. A histogram 
of the 754 Sh measurements for this star is plotted in Figure 8. The FWHM is 0.008 and 
represents the empirical uncertainty in our single measurement precision for Sh at Lick 
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Observatory. This single measurment precision is about five times larger than the Keck Shk 
values, but the median Sh value from Lick is 0.1627 ± 0.008 and agrees with the Keck value 
Shk = 0.167 ± 0.002 within uncertainties. 



Results 



In total, more than 44,000 spectra of 2630 stars were analyzed and calibrated to the 
Mt. Wilson Shk activity index. The Shk and B — V values were used to calcul ate log.R HK , 



a me asure of the flux in the Ca II lines relative to the basal photospheric value (jNoyes et al. 



1984 ). While Shk has a functional dependence on spectral type, logi? HK removes the basal 
component of core emission and only varies with chromospheric activity. The \o&R! WK pa- 



rame t er has been calibrated to rotat ion periods and stellar ages for FGK stars (jNoyes et al. 
19841 ; iMamaiek fc Hillenbrand! boost ). 



The main disadvantage of logi? HK is that stars with B — V less than 0.4 or greater 
than 1.0 were not included in the original calibr ations because main sequence stars outside 
this r ange were not observable in distant clusters (jNoyes et al.lll984l ; iMamajek &: Hillenbrand 
20081 ). The derived quantities logi? HK , -^rot, and stellar age are most secure for main sequence 
stars that have log.R HK between -4.0 and -5.1 and B — V values between 0.4 and 1.0. 

Table [1] contains a summary of star names, B — Vcolor, the median values of Shk and 
logi? HK height (8My) above the main sequence, P rot (days), Age (Gyr) and the Observatory 
where the measurements were obtained (Keck or Lick). A stub of this table is provided in 
the printed version of this paper and the complete table is available in the online version of 
this paper. 

For evolved stars and stars outside the color range 0.4 < B — V < 1 .0, the derived param- 
eters l og R' RK P ro t and Age were not calibrated by lNoyes et al.l (119841 ) or lMamajek fc Hillenbrand 
(120081 ) and are not included here. The complete time series data from Keck and Lick Observa- 
tories are listed in Table [2J The spectra from Keck Observatory were collected between 2004 



Augus t and 2010 August. These values can be appended to results contained in lWright et al. 



(j2004j ) so that for many stars, the measurements approach the span of an entire activity cycle 
(e.g., an eleven year solar cycle). The time series observations in Table [2] (online data only) 
from Lick Observatory span 2002 through 2009 December. 
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4.1. Outliers 



A few specific stars have va lues very different f r om th e published values of lDuncan et al. 



f |199lh and IWright et al.l (1200 



a typographical error (IWright 



In IWright et al.l (120041 ) . the S HK value for HD 38392 has 



2009|), HD 104958 has an S HK of 0.07 which is likely a typo- 



graphical error (bec ause it is non-physica l) and HD 220339 has an Shk value of 0.000, which 
is also non- physical IWright et al.l (120041 ) . We also note that HD 195019 is a binary with 



S-values for the A and B components. 



4.2. Chromospheric activity of main sequence stars 

Stars in the CPS sample are defined here as main sequence stars if the absolute visual 
magnitude is within 1.5 magnitudes of the main sequence. The simultaneous measurements 
of activity and precision velocities are an ideal data set to search for correlations between 
activity and spurious radial velocity jitter. Unfortunately, log 5 HK nas n °t been calibrated 
for evolved or late type stars. An activity metric is desirable because these stars are being 
monitored for exoplanet surveys and we would like to understand the impact of chromo- 
spheric activity on radial velocities for evolved and late-type stars (as well as main sequence 
stars) to quickly assess whether velocity scatter is likely to be astrophysical or dynamical in 
origin. 

As a proxy for log R' KK , which accounts for the basal (rotation independent) photospheric 
flux near the H & K lines, we define baseline activity values, Sbl, by fitting the floor in Shk 
for main sequence stars as a function of 5 — V color: 



Sbl =2.7 - 16.19(5 - V) + 36.22(5 - V) 2 - 27.54(5 - Vf - 14.39(5 - V) 4 (5) 
+ 34.97(5 - Vf - 18.71(5 - Vf + 3.17(5 - V) 7 

The polynomial fit described by Eqn 5 for Sbl is shown in Figure 9 as a dashed red 
line and was obtained by selecting the tenth percentile (i.e., inactive) Shk values in twelve 
5 — V bins (0.1 magnitude in width). Statistical error bars were calculated by dividing the 
tenth percentile activity value by the square root of the number of points in each bin. The 
error bars were used to weight the polynomial fit to the baseline Shk- Because of the large 
stellar sample, we were able to carry out this empirical fit, rather than relying on theoretical 
estimates to remove the basal emission in the Ca II line core. It is reasonable to assume that 
each of the color bins contains a fraction of low activity stars. 
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The Sbl values for main sequence stars begin to rise redward of B — V = 1.0, expected 
in part because of a decrease in continuum flux for redder stars. Figure 10 demonstrates the 
genuine range in Ca II core emission, plotting Ca II lines for three stars with B — V ~ 1.5 
but different Shk values. The lower env elope for Shk turns over at B — V of 1.4. This is 
consistent with activity measurements of iRauscher fc Marcyl (120061 ) who measured line core 
equivalent widths and show that the luminosity in the Ca II lines decreases by a fac tor of 
~ 3 with decreasing mass from K7 to M5V. Likewise, iGizis. Reid. fc Hawleyi (120021 ) mea- 
sured Ha EW as an indicator of chromospheric activity in M dwarfs and found that the 
E W of absorption line s decre ases with increasing spectral type from K7V to M4V. Models 
of ICram &: Giampapal (119871 ) show that M dwarfs with the strongest Ha absorption have 
moder ately active chromospheres , while weak Ha correlates with weak chromospheric ac- 
tivity. ICram fc Giampapal ( 119871 ) further show that activity in M dwarfs may be relatively 
constant over time, in contrast with solar type stars with strong activity cycles. 



4.3. Chromospheric activity in subgiants 



The activity measurements in this pa per include a substantial su bset of 234 subgiants, 
included a su rvey that targets su bgiants (IJohnson et al.l 120071 . 120101 ) and contained in the 
N2K survey (IFischer et al.l 120051 ) . The Shk values for these stars, plotted in Figure 11, 
show a modest delcine in the median Shk value as a function of B — V. A typical Shk 
value for subgiants is 0.1 < Shk < 0.2. However, 10% of the 234 stars show moderate to 
strong chromospheric activity, with Shk values greater tha n 0.3. We have identifie d some 
of t hese active stars as sp e ctrosc opic binaries (blue dots). iGunn. Mitrou fc Doyld (119981 ) 
and IGizis. Reid. fc Hawleyi (120021 ) have noted a correlation between chromospheric activity 
and orbital periods for both main sequence and evolved stars. In Figure 12, the median 
Shk value is plotted as a function o f vsini. Stars w ith higher v sinz are also more active, 
consistent with the idea proposed by iGillilandl (119861 ) that activity may be a transient phase 
for subgiants. However, the sprinkling of ac tive stars is found acr oss the subgiant branch, 
not just at the G0IV boundary identified by lSimon fc Drake! (119891 ). 



5. Activity and jitter 

The primary motivation for studying chromospheric emission in planet search stars is 
to monitor potential sources of astrophysical noise that can impact the radial velocities. 
We calculate the difference between the median Shk value and Sbl as a AS for each star. 
This AS is an empirical parameterization for the excess activity of stars at all B — V colors, 
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including late K and early M dwarfs. Because the photon-limited velocity precision at Keck 
is better than Lick, we restricted our assessment of velocity jitter to stars observed with 
HIRES at Keck. 

We divided the target stars into four B — V groups and plotted velocity rms as a function 
of AS* in Figure 13. We did not remove any linear trends or periodic velocities, even from 
stars with known planets, so some of the scatter above the lower envelope of velocity rms 
is due to dynamical velocities. We define jitter as the quadrature difference of velocity rms 
minus the formal internal errors (Eqn 6). We emphasize that since velocity scatter includes 
instrumental and analysis errors, this definition of jitter includes systematic instrumental 
errors as well as astrophysical noise. Furthermore, the (astrophysical or instrumental) jitter 
for any particular star may be higher than the tenth percentile velocity rms floor. 

Jitter = V rms 2 — a 2 (6) 

Figure 13 (top, left) contains data for 259 late F to mid G-type main sequence stars 
with 0.4 < B — V < 0.7. The minimum jitter as a function of AS is fit to the lower tenth 
percentile in velocity rms with a line that is plotted in red in the top left panel of Figure 
13. For chromospherically quiet stars, the level of jitter begins at 2.3 m s~ 1 and increases 
rapidly as a function of AS (Eqn 7). Note that the AS range is smaller for the blue stars 
because the photospheres are brighter at blue wavelengths and there is less contrast in the 
Ca II lines. The stars on the CPS sample in this bluest B — V color bin exhibit weaker 
Ca II core emission, however, the log R' HK value for these stars climbs steeply as a function 
of increasing Shk- In our experience, these stars have the strongest velocity jitter of the four 
groups considered. 

The second group consists of 218 main sequence stars in the color range 0.7 < B — V < 
1.0. The tenth percentile floor in velocity rms is fit by Eqn 8 and plotted as a red line in the 
top right panel of Figure 13. The minimum jitter for chromospherically quiet stars in this 
B — V color range is 2.1 m s -1 and jitter increases with increasing activity, AS. 

The third color bin, 1.0 < B — V < 1.3 contains 118 stars. These stars have a minimum 
noise floor of 1.6 m s _1 and Eqn 9 does not have a dependence on AS. Even when the 
emission in the core of the Ca II lines is significant, it is possible to find stars in this color 
bin that have low jitter. 

The noise floor for the 89 reddest stars (bottom right panel of Figure 13) begins at 2.14 
m s -1 and Eqn 10 shows that jitter increases slowly with AS . It is possible that the noise 
for this group of stars is affected by analysis errors, since the deconvolution is particularly 
challenging for spectra of late type stars. 
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Jitter = 2.3 + 17.4 * AS m s' 1 (OA < B-V < 0.7) (7) 

Jitter = 2.1 + 4.7* AS m s" 1 (0.7 < B-V < 1.0) (8) 

Jitter = 1.6 - 0.003* AS m s" 1 (1.0 < B-V < 1.3) (9) 

Jitter = 2.1 + 2.7 * AS m s" 1 (1.3 < B-V < 1.6) (10) 

We have also carried out an analysis of activity-correlated jitter for subgiant stars ob- 
served at Keck, defined here as stars that have a 5My of at least 1.5 magnitudes above 
the main sequence. Eqn 11 provides our fit to the bottom tenth percentile activity values. 
Ninety percent of subgiants have low chromospheric activity, so the fit to the floor of the 
velocity rms in Figure 14 spans a relatively narrow range of AS. The jitter for subgiants is 
defined in Eqn 12. 

Sbl = 0.2 — 0.07 * B — V (11) 
Jitter = 4.2 + 3.8 * AS m s" 1 (5M V > 1.5) (12) 



6. Discussion 

The Ca II lines are good indicators of chromospheric activity for main sequence stars. 
Here, we present measurements of this emission, parameterized as Shk values, for ~ 2600 
stars observed at Keck Observatory and Lick Observatory. The time series of activity mea- 
surements at Keck Observatory date back to 2004 August and the time series data for Lick 
began in 2002 when the spectral format was extended to include the Ca II lines. 

To put Shk on a continuous time baseline, we calibrated our measurement s to the values 



obtained from the long-standing Mt Wilson Observatory H & K program (IWilsonl 11968 



Duncan et al.lll99ll ) using a set of stars common to both projects. Although we developed 



a new approach for diffe rentially measuring Shk, these time series activity measurements 



may be appended to the IWright et al.l (120041 ) values, which were also calibrated to the Mt 



Wilson H&K program. We find an rms of 11% in the difference between our Shk values 
and those on the long term Mt Wilson H&K project. Since the rms of Shk values for 
chromospherically quiet stars is only 1% for Keck stars and 5% for Lick stars, we expect 
that long term chromospheric variability contributes to the 11% rms in our calibration. Our 
measurement of Shk includes a significant sample of stars redward of B — V = 1.0. 

The Mt Wilson H&K project surveys the brighter F, G and K-type main sequence 
stars and subsequent calibrations to logi? HK , P rot and ages were restricted to these spectral 
types. With the benefit of a large statistical sample of stars, we define the lower envelope 
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of activity, Sbl, as a function of B — V for main sequence and subgiant stars over the color 
range 0.4 < B — V < 1.6. We define AS as the difference between the median 5hk f° r 
each star and Sbl- We then evaluated velocity jitter for our stars by subtracting the mean 
internal error from the rms of the velocities for each star. 

Velocity "jitter" therefore characterizes the extent to which the velocity rms deviates 
from our la internal errors. Jitter includes instrumental and analysis errors as well as 
photon shot noise, unresolved dynamical velocity shifts, and photospheric noise, and it is 
not always clear which of these is the dominant term. In programs with uniform analysis 
precision where a constant signal-to-noise is obtained for all observations (the case for CPS 
observations of bright stars at Keck), it is possible to treat jitter as a free parameter in 
the Keplerian model. However, both Frequentist and Bayesian models can benefit from 
physical priors since additional unresolved low amplitude planet signals will inflate residuals 
to the Keplerian fit. It is likely that jitter will be over-estimated (and more ambiguous) if 
velocities are obtained with a spectrometer where hardware components (e.g., the detector 
and controller) have changed, or if data sets have been obtained with more than one telescope, 
or if the signal-to-noise changes significantly from observation to observation. 

A key conclusion from this work is that K dwarfs in the color range 1.0 < B — V < 
1.3 have the lowest level of velocity jitter. Importantly, jitter appears to be decoupled 
from chromospheric activity in these stars; it is a constant 1.6 m s _1 independent of the 
strength of emission in the Ca II line cores. The lack of correlation between observed jitter 
and chromospheric activity in K dwarfs suggests that jitter for these stars is dominated by 
instrumental or analysis errors and not astrophysical noise sources. This demonstrates two 
important points: (1) on all timescales since 2004, the systematic instrumental errors in CPS 
observations at Keck are no more than 1.6 m s _1 and (2) the mid to late K dwarfs represent 
a sweet spot for optimal detectability in Doppler exoplanet searches. The astrophysical 
contribution to velocity noise in K dwarfs could be significantly less than 1 m s _1 . 

Relative to Kdwarfs, we also show that jitter increases as a function oi B — V for chromo- 
spherically inactive stars. The lowest activity stars in the bluer color bin of 0.7 < B — V < 1.0 
and in the redder color bin of 1.3 < B — V < 1.6 exhibit a measurably higher minimum jit- 
ter of about 2.1 m s _1 and a weak dependence on activity (AS). The minimum jitter 
increases slightly to 2.3 m s _1 for chromospherically inactive late F and early G dwarfs 
(0.4 < B — V < 0.7) and there is a strong sensitivity to AS. 

If we assume that most of the observed K dwarf jitter is instrumental in origin (i.e., 
we assume that the astrophysical contribution to jitter is essentially zero for K dwarfs), 
then stars in other color bins exhibit a quadrature sum of 1.6 m s _1 instrumental noise plus 
additional intrinsic astrophysical noise. Thus, excluding instrumental contributions, the 
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astrophysical noise floor imposed by chromospherically inactive late F and early G dwarfs 
appears to be at least a/2.3 2 — 1.6 2 = 1.7 m s -1 and we expect that it will be difficult, even 
with the best Doppler precision, to beat down this noise floor. It may be possible to reduce 
the impact of stellar jitter in late F and G-type dwarfs by confining Doppler searches to 
stars in the bottom quartile of activity level and by observing the stars very frequently (i.e., 
with a high cadence strategy). However, in general, solar twins are not optimal targets for 
Doppler exoplanet searches that aim to detect very low mass planets. At best, sunlike stars 
will certainly require more data points to average down a noise floor that starts out higher 
than the astrophysical noise from K dwarfs. 

By this same reasoning, mid- to late-G dwarfs with 0.7 < B — V < 1.0 impose an 
astrophysical noise floor that is ~ 1.3 m s^ 1 above that of K dwarfs. It is possible that 
later type stars with 1.3 < B — V < 1.6 may have greater jitter because of challenges with 
deconvolution (i.e., a noise source that is not astrophysical). A dependence of rms velocities 
on AS would argue for astrophysical noise, however the functional dependence is quite weak 
and the number of stars with larger AS puts us in the realm of small number statistics. 

Since there is a correlation between AS" and the rms velocity floor in stars, one might 
hope that variations in individual observation s of 5hk would have a nearly one-to-one cor- 



relation with Doppler velocity measurements. iQueloz et al.l ( 1200 ll ) reported the remarkable 



case of HD 166435, a young GOV star with quasi-period radial velocities that were correlated 
to line bisector variations and Ca II emission. In this unusual case, the star seemed to have 
persistant active longitudes where spots were regenerated for almost two years, producing 
a coherent signature of magnetic activity. However, stars do not generally show such clear 
correlations between activity and velocity. 

A more typical example is HD 143714, a metal-rich GOV star from the N2K program 
with B — V = 0.61 and a median S'hk = 0.24. According to Eqn 5, Sbl = 0.147 for a 
star of this color, so AS = 0.093. Using Eqn 7, the expected jitter for HD 143714 is 3.9 
m s -1 . By the end of 2008, the velocities for HD 143714 showed a significant rms of almost 
10 m s _1 and were periodic (Figure 16). However, we observed an apparent correlation 
between the Doppler velocities and activity (Figure 17). When the velocities were detrended 
using the linear fit to activity, the periodicity disappeared and the rms dropped to 6.2 m s _1 . 
We apparently observed an epoch where chromospheric activity and radial velocities were 
correlated and the level of velocity rms exceeded the predicted jitter level. However, the 
activity correlation disappeared. Velocity measurements obtained after 2009 January are no 
longer correlated with activity and have an rms of 4.9 m s _1 closer to the predicted jitter. 
Similar correlations between 5hk and velocities have been observed in only a few other CPS 
stars. Like HD 166435 and HD 143174, the stars are late F or early G dwarfs and the 
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correlations come and go. The distribution of jitter values for CPS stars is summarized in 
the HR diagram in Figure 15. 

Subgiants generally exhibit lower chromospheric activity than main sequence stars, 
however these stars apparently have a jitter floor starting at about 4 m s _1 , significantly 
higher than the jitter level of main sequence stars with similar B — V. As the stars expand, 
the depth of the convective zone increases, braking rotation and decreasing the convective 
turnover timescale. Consistent with this, ninety percent of the subgiants analyzed here 
have low chromospheric activity. However, ten percent of the subgiants have moderate to 
strong chromospheric activity that is correlated with rotational velocity (v sin i) or the pres- 
ence of a clo se binary c ompa nion. These active subgiants are reminicent of the F8IV star, 
HD 125840. iGillilandl ( 119861 ) explained the unusual activity of this star as the onset of a 
convective dynamo d uring the evolution of t he m assive progenitor star. Activity in subgiants 
has been observed by iGray fc Nagarl (119851 ) and ISimon fc Drake! (119891 ) for subgiants hotter 
than GOIV. However, the active subgiants here are not limited to this boundary region; they 
are found across the subgiant branch. This range of subgiant activity may reflect a range of 
conditions in the progenitor star, such as stellar mass or rotation, that transitions through 
a convective dynamo stage at different points on the subgiant branch. Beca use only a small 
fractio n of subgiants are active, this transition phase is likely short-lived. ISimon &: Drake 
( 119891 ) claim that activity in subgiants is restricted to 100 Myr; the time for early F-type 
stars to evolve from the main sequence to a GOIV star, implying that the convective zones 
of most subgiants does not sustain convective dynamos for long. Although the phenomenon 
of activity in subgiants is uncommon, the existence of active subgiants with spectral types 
later than GOIV sugge sts that this phase can persist longer than the 100 Myr estimate by 
Simon fe Drakel f)l989h . 



We thank Jason T. Wright for useful conversations. We gratefully acknowledge years of 
work on the Keck planet search program by G. W. Marcy, R. P. Butler and S. Vogt. We also 
acknowlege the many observers who collected observations and contributed to the Doppler 
analysis: J. A. Johnson, A. W. Howard, Julien Spronck, Kelsey Clubb, John Brewer, K.Peek, 
J. Anderson, J. A. Valenti, N. Piskunov. We gratefully acknowledge the dedication and 
support of the Keck Observatory staff, in particular Grant Hill and Scott Dahm for support 
with HIRES and Greg Wirth for supporting remote observations. D.A.F. acknowledges 
research support from NASA grant NNX08AF42G and NSF AST 1036283. We thank the 
NASA Exoplanet Science Institute (NExScI) for support through the KPDA program. We 
thank the NASA and NOAO Telescope assignment committees for allocations of telescope 
time. The authors extend thanks to those of Hawaiian ancestry on whose sacred mountain 
of Mauna Kea we are privileged to be guests. Without their kind hospitality, the Keck 



- 16 - 



observations presented here would not have been possible. This research has made use of 
the SIMBAD database, operated at CDS, Strasbourg, France, and of NASA's Astrophysics 
Data System Bibliographic Services. 



-17- 



REFERENCES 

Apps, K. et al. 2010, PASP, in press 
Baliunas, S. L., et al. 1995, ApJ, 438, 269 

Bonfils, X., Mayor, M., Delfosse, X., Forveille, T., Gillon, M., Perrier, C, Udry, S., Bouchy, 
F., Lovis, C, Pepe, F., Queloz, D., Santos, N.C., Bertaux, J.-L. 2007, A&A, 474, 293 

Cram, L.E., Giampapa, M.S. 1987, ApJ, 323, 316 

Duncan, D. K., et al. 1991, ApJS, 76, 383 

Fischer, D. A., et al.2005, ApJ, 633, 465 

Gilliland, R. 1986, Lecture Notes in Physics, Springer Berline/Heidelberg 254, 177-180. 
Gizis, J.E, Reid, I.N., Hawley, S.L. 2002, AJ, 123, 3356 
Gray, D., Nagar, ApJ298, 756 

Gunn, A. G., Mitrou, C. K., Doyle, J. G. 1998, MNRAS296, 150 

Hall, J., Lockwood, G. W. & Skiff, Brian A. 2007, A J, 133, 862H 

Henry, T. J., Soderblom, D. R., Donahue, R. A., & Baliunas, S. L. 1996, AJ, 111, 439 

Henry, G. , Donahue, R. , & Baliunas, S. 2002, ApJ, 577, Llll 

Jenkins, J., Jones, H.R.A., Tinney, C. G., Butler, R.P., McCarthy, C, Marcy, G.W., Pinfield, 
D. J., Carter, B.D., Penny, A.J. 2006, MNRAS, 372, 163 

Johnson, J. A., Marcy, G.W., Fischer, D.A., Laughlin,G., Butler, R.P, Henry, G.W., Valenti, 
J.A., Ford, E. B., Vogt, S.S. & Wright, J.T. 2006, ApJ, 647, 600 

Johnson, J. A., Fischer, D.A., Marcy, G.W., Wright, J. T., Driscoll, P., Butler, R. P., Hekker, 
S. Reffert, S., Vogt, S. S. 2007, ApJ, 665, 785 

Johnson, J.A., Howard, A.W., Bowler, B.P, Henry, G.W., Marcy, G.W., Wright, J. T., 
Fischer, D.A., & Isaacson, H. 2010, PASP, 122, 701 

Kraft, R. 1967, ApJ150, 551 

Mamajek, E. E.; Hillenbrand, L. A. 2008, ApJ, 687,1264 



Marcy, G. W., Butler, R. P., Fischer, D. A., Vogt, S. S., Wright, J. T., Tinney, C. G., Jones, 
H. R. A. 2005 PThPS, 158, 24 

Middelkoop, F. 1982, Al07, 31 

Noyes, R. W., Hartmann, L., Baliunas, S. L., Duncan, D. K. & Vaughan, A. H. 1984, ApJ, 
279, 763 

Paulson, D. B., Cochran, W. D., Hatzes, A. P. 2004, AJ, 127, 3579 

Queloz, D., Henry, G. W., Sivan, J. P., Baliunas, S. L., Beuzit, J. L., Donahue, R. A., Mayor, 
M., Naef, D., Perrier, C. & Udry, S. 2001, A&A, 379, 279 

Rauscher, E., & Marcy, G.W. 2006, PASP, 118, 617 

Saar, S. H., Butler, R. P. & Marcy, G. W. 1998, ApJ, 498, 153 

Saar, S. H. & Fischer, D. A. 2000, ApJ, 534, 105 

Spronck, J., Fischer, D. A. 2010 (in prep) 

Shine, R. & Linsky J. 1974, Solar Physics, 39, 49S 

Simon, T., Herbig, G. H., & Boesgaard, A. M. 1985, ApJ293, 551 

Simon, T. Drake, S. A., 1989, ApJ346, 303 

Skumanich, A. 1972, ApJ171, 565 

Soderblom, D. 1983, ApJS53, 1 

Strassmeier, K., Washuettl, A., Granzer, T., Scheck, M., & Weber, M. 2000,A&AS, 142, 275 

Valenti, J. A., & Fischer, D.A. 2005, ApJS, 159, 141 

Vogt, S. S. et al. 1994, SPIE, 2198, 362. 

Vogt, S. S. 1987, PASP, 99, 1214 

Wilson, O. C. 1963, ApJ, 138, 832 

Wilson, O. C 1968, ApJ, 153, 221 

Wright, J. T. (private communication) 

Wright, J T. et al. 2004, ApJ, 128, 1273 (erratum, 129, 1776) 



- 19 - 



Wright, J. T. 2005, PASP, 117, 657 



This preprint was prepared with the A AS IATgX macros v5.2. 



-20 - 




3890 3900 3910 3920 3930 3940 3950 
Wavelength (Angstroms) 



Fig. 1. — The echelle order from Keck/HIRES containing the Ca II K line for the G8V star, 
r Ceti. Emission in the core of the Ca II K line is weighted by the indicated normalized 
triangle, giving greatest weight to flux in the core of the line. The Ca II flux is measured 
relative to continuum flux in the indicated V band. 
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Fig. 2. — The echelle order from Keck HIRES with the Ca II H line near the center of 
the echelle order. The K line toward the edge of the detector where the blaze function is 
dropping was not used because of the low SNR in the line and continuum flux. The triangular 
weighting function is indicated on the plot for the H line. The wings of the H & K lines in 
this order suppress the continuum so that the R continuum band is defined in an adjacent 
order. 
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Fig. 3.— The 20 A R band was denned in an adjacent order to monitor photospheric 
continuum redward of the Ca II H & K lines. This spectrum in this figure is a Keck/HIRES 
spectrum of r Ceti. 



-23 - 




Fig. 4. — The calibration of co efficients for Sw 



was made using 151 stars observed in common 



with the Mt. Wilson project ( iDuncan et al.lll99ll ). The free parameters (coefficients) were 



determined with a Levenberg Marquardt algorithm that minimized the difference between 
Keck and Mt. Wilson Shk values. The calibration yielded an RMS scatter of 11% for 
stars of all spectral types and all activity levels. Variability in chromospheric activity also 
contributes to some of the observed RMS scatter. 
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Fig. 5. — A histogram of Shk for the chromospherically inactive star, r Ceti. The median 
S'hk value from Keck is 0.167 and the FWHM of the distribution is roughly 0.009 or 1%. 
This provides an empirical estimate for the errors assigned to activity measurements. 
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Fig. 6. — Only one echelle order contains adequate signal for analysis of Shk at Lick Ob- 
servatory. The "continuum" section as well as the weighting function for the H line core is 
indicated in this spectrum of r Ceti. 
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Duncan (1991) S-values 



Fig. 7. — The Sr values measured at Lick are based only on emission in the Ca II H line. 
There were 83 stars in common between the Lick program and the Mt. Wilson H & K 
survey that were used for the calibration. An RMS scatter of 11% was found when fitting 
the coefficients for Sr at Lick. 



0.14 0.15 0.16 0.17 0.18 0.19 

Lick S-values for Tau Ceti 

Fig. 8. — The histogram of Su values for r Ceti from Lick Observatory has a median Su 
value of 0.163. The FWHM of the measurement distribution is roughly 0.008, suggesting a 
single measurement precision of about 5% for Su. 
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Fig. 9. — The 5hk values for main sequence stars are plotted as function of B — V color. The 
lower envelope of the Shk values is defined as the tenth percentile in activity (red dashed 
line) and adopted as a basal activity level, Sbl- There is a decrease in Sbl for stars redward 
of B — V = 1.4 that reflects a drop in chromospheric activity. 
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Fig. 10. — Emission lines are shown for three stars with approximately the same B — V and 
different activity levels that demonstrate a range in chromospheric activity. 
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Fig. 11. — The Shk values of the vast majority of subgiants is consistent with low chromo- 
spheric activity. About ten percent of subgiants are active. Some of the active subgiants are 
rapid rotators or close binaries. 
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Fig. 12. — The chromospheric activity in subgiants is plotted as a function of v sin i. Rotation 
and activity appear to be correlated. The presence of a close binary stellar companion is 
also correlated with chromospheric activity. 
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Fig. 13. — Radial velocities are plotted for main sequence stars as a function of Sbl- The 
left plot shows stars with 0.4 < B — V < 1.0 and the right plot shows data for stars with 
1.0 < B — V < 1.6. The floor of the velocity rms is adopted as jitter for these stars (which 
may include instrumental systematic errors). Overall, the redder stars exhibit a lower level 
of jitter than blue stars. High Ca II emission stars are rare among the bluer stars, but the 
velocity jitter is steeper for these stars. 
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Fig. 14. — Radial velocities are plotted for the IJohnson et al.l ( 120101 ) sample of subgiants. 
Most subgiants are chromospherically inactive, so it is difficult to quantify the slope of jitter 
for these stars. However the subgiants appear to be only modestly more active than main 
sequence stars. 
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Fig. 15. — The HR diagram for CPS stars. The symbol color is scaled according to the jitter 
legend and symbol size is scaled to v sini. The lower main sequence stars have modest jitter, 
while FGK stars show a broader range in background noise. Most stars on the subgiant 
branch have low activity and low intrinsic jitter, however about ten percent of subgiants are 
chromospherically active. 
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Fig. 16. — Radial velocities for HD 143174 exhibited velocity rms of almost 10 m s -1 and clear 
structure in the time series Doppler measurements that was initally well-fit by a Keplerian 
model. The activity level for this B — V = 0.61 star has a AS = 0.093 and expected jitter 
of just 3.9 m s _1 . 
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Fig. 17. — HD 143174 shows a suspicious correlation between chromospheric activity index, 
Shk and radial velocities. 
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Table 1. S- values and Derived Quantities 



Star 






Median 


Jitter 


Median 




Prot 






name 3 


HIP 


B-V 


S 


TTIS _ 1 


log Bin, 
1KJ b - 1 ^HK 


AM, 


( davs) 




Obsv 


GJ 26 




1.540 


0.895 


2.844 




0.17 






K 


38A 


473 


1.410 


1.775 


4.542 




1.01 






K 


38B 




1.410 


1.747 


4.475 




0.15 






K 


GJ 47 




1.390 


1.074 


2.885 




-1.4 






K 


GJ 48 


4856 


1.460 


0.907 


2.483 




-0.9 






K 


GJ 49 


4872 


1.500 


2.947 


7.510 




0.51 






K 


GJ 83.1 




1.820 


14.72 


0.000 




0.76 






K 


GJ 87 


10279 


1.430 


0.634 


2.100 




-0.8 






K 


105 


490 


0.600 


0.380 


6.581 


-4.344 


-0.0 


3 


0.19 


K 


GJ 105B 




1.620 


0.843 


3.717 




-0.1 






K 


GJ 107B 




0.930 


1.766 


9.987 




-3.0 






K 


GJ 109 


12781 


1.550 


0.814 


2.736 




-0.2 






K 


166 


544 


0.750 


0.434 


3.586 


-4.384 


0.00 


7 


0.26 


K 



a Star names are HD unless otherwise given. Values of B-V are from the Hipparcos Catalog. 
Derived quantities are only given for stars in the range —4.0 < logi?^ K < —5.1. Obsv is K 
for Keck Observatory or L for Lick Observatory. 
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Table 2. S-values and Julian Date - 2440000 



Star 


JD 










name HIP 


- 2440000. 


S value 




SNR 


Obsv 


GJ 26 


13238.978 


0.686 




17 


K 


GJ 26 


13338.806 


0.614 




13 


K 


GJ 26 


13369.738 


0.661 




15 


K 


GJ 26 


13426.728 


0.598 




8 


K 


GJ 26 


13723.735 


1.045 




9 


K 


GJ 26 


13961.013 


0.939 




14 


K 


GJ 26 


13981.934 


0.939 




15 


K 


GJ 26 


13982.964 


0.845 




14 


K 



